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1. Introduction 

The theory of noise in quantum transport in mesoscopic systems is a very active field of 
research IH [2]- In addition to the first few moments of the transmitted charge the full 
probability distribution can be studied, called full counting statistics (FCS). This was 
first done in a publication by Levitov and Lesovik [3] where noninteracting fermions were 
treated using the cumulant generating function. The system usually studied consists 
of a finite "dot" -region connected to M leads which initially are separated from the 
dot region and have different chemical potentials [21 HI [5|, [6]. After connecting the 
subsystems the time evolution of the particle transfer between the leads is studied. In 
order to avoid mathematical subleties it is useful to start with a finite total number of 
particles A^tot which can be achieved using leads of finite extent. The thermodynamic 
limit is performed only at a later stage. For the lattice models 

with a finite number of states at each lattice site studied in this paper this implies 
also a finite number Nh of the dimension of the Hilbert space of a single fermion. 

For an initial state which is a Slater determinant with iVtot fermions the 
characteristic function g{t) for noninteracting fermions is a Atot x Atot determinant. 
After averaging over a grand canonical ensemble an expression for g{t) in terms of a 
Nff X Nh determinant can be derived [H [3 [8] . In some publications this result is called 
"Levitov-Lesovik formula" 0. This expression is the formal starting point for the actual 
calculation of the characteristic function. It consists of two steps. The first one is to 
calculate the time dependence of the one-particle projection operators Pa(t) onto lead 
a. For finite times exact results can be obtained numerically [8J. In the long time 
limit an accurate analytical approximation can be given in terms of scattering states 
after performing the thermodynamic limit. In the second step the determinant over 
the one-particle Hilbert space has to be calculated. After the thermodynamic limit 
has been performed this is an infinite dimensional determinant and mathematical care 
is necessary [10]. Various approaches have been proposed for the evaluation of the 
determinant. Muzykantskii and Adamov [Tl] used methods from the theory of singular 
integral equations to proceed for M = 2. In the long time limit they obtained the 
leading term for lng{t) (linear in t) by the exact solution of a matrix Riemann-Hilbert 
problem. Their approach provided the first explicit derivation of the result presented 
by Levitov and Lesovik [3]. Alternatively one can use a formal power series expansion 
of the logarithm of the determinant [H]- The term linear in time can then easily be 
identified and the infinite dimensional determinant can be reduced to an energy integral 
over the logarithm of a M x M determinant. A third approach used Szego's theorem 
from the theory of Toeplitz matrices to obtain the term linear in time [12] . 

For two leads (M = 2) the contribution to ln(7(t) linear in t vanishes in the case 
of perfect transmission and subleading terms logarithmic in time have to be considered. 
This was studied analytically by Muzykantskii and Adamov [11] by an approximate 
solution of a more complicated Riemann-Hilbert problem as well as by extensions of 
Szego's theorem [12]. Numerical as well as analytical results in agreement with these 
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findings were presented by one of us|8j. 

In tliis paper we generalize and simplify the derivation using the formal power series 
expansion for In g{t). After obtaining the general linear in t contribution the M x M 
determinant for the joint probability distribution is examined in detail. The M x M 
matrix of which one has to evaluate the determinant is written as sum of the unit matrix 
and a second matrix. This allows to read off the general expression for the characteristic 
function for the joint probability distribution of two or more observed charged transfers 
for arbitrary values of M. 

Applications to Y-junctions (M = 3) and dots with M = 4 leads are discussed. 

2. General formulation 

In the following we consider a system which consists of a finite "dot" -region described 
by the Hamiltonian H^"^ and M leads with the Hamiltonians Ho a with a = 1,..,M. 
The leads are initially separated from the dot region. The number of electrons in the 
initial state are iVg"* and No^a- We assume the initial state |$(0)) to be an eigenstate 
of H^"^ and the i/o,a 

|<|.(0)) = |Ef°") ® ® ... ® \E^"^-'). (1) 

The time evolution for times greater than zero is described by the Hamiltonian 

H = + 5^//o,a + J^K = i^o + ^ • (2) 

a a 

The term V which couples the leads with the dot region will be specified later. The 
probability distribution that Qi electrons are transfered to lead 1, Q2 electrons are 
transferred to lead 2 etc. after time t when the subsystems are separated again, is given 

by 

M 

w{t, {Q}) = imiU^iQa - {K - No,a)]\m) (3) 

a=l 

= jiy^JdX^--->^Me~'^^'^^'9it,{X}). 

Here A/l ist the particle number operator of the lead a and g{t, {A}) is the characteristic 
function. With the particle number operators Ma{t) in the Heisenberg picture g is given 
by 

g{t^ {A}) = (<l>(0)|e'^^"^"We-'^^"^''|<l>(0)) . (4) 

The fact that the initial state is assumed to be an eigenstate of the particle number 
operators was used. This expression can be easily generalized to initial statistical 
operators of the type po = Po°* ® Po^^ ® ... (8> Pq ^^ ■ An important example are initial 
grand canonical subensembles with different temperatures and chemical potentials 
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and Pq°* of the same type. Then po has the generahzed canonical form pq = e ^°/Zq. 
Averaging yields for the characteristic function 
g{t, {A}) = (e*EAaA^„(^)g-^EA.Ar.^ ^ 

where (...) denotes the averaging with the statistical operator po- This result is also 
valid for interacting fermions. 

For noninteracting fermions the expectation value can be simplified using Klich's 
trace formula (TJ [H] 

Tr(e^e^) =det(l + e'^e') , (7) 

where A and B are arbitrary one particle operators in Fock space and a and b are 
the corresponding operators in the Hilbert space of a single particle. Therefore the 
characteristic function can be expressed as a determinant in the one particle Hilbert 
space 

g{t, {A}) = det [1 + [e'^a^aPa{t)^-ij:^XaPa _ ^^j 

^det[l + 6(t,{A})] , (8) 

where Pa is the projection operator in the Hilbert space of one particle on the states of 
the a-th lead and uq = {e^° + l)~^is the Fermi operator. It is determined by the Fermi 
functions describing the initial state. In order to obtain joint probability distributions 
for arbitrary times the first step is the (numerical) calculation of the Heisenberg operator 

Pa{t). 

In order to study the long time behaviour it is useful to introduce the current 
operators ja = [Pa,h]/i, where h is the Hamiltonian in the Hilbert space of a single 
particle, and write Pa(t) as [8] 

Pa{t) = Pa{0) + fja{t')dt' = + 5P,(t) . (9) 

The operator b in Eq. (IHl) can then be expressed as 

b{t, {A}) = ^(e'^" - l)5Pa(t)e-*^-'^"'-^'"'no . (10) 

a 

In order to avoid reflections from the ends of the leads (far away from the dot) in 
the long time limit, the thermodynamic limit has to be taken first. In this limit the 
discrete energy spectrum of the initially occupied standing wave states \ej,a) becomes 
continuous and the trace in the one-particle Hilbert space involves an energy integration 
for the lead states 

tr6 = 5^(6f>|en + E / de{e,a\b\e,a) , (11) 

i a 

with the normalization (e, a|e',a') = Saa'S{e — e'). Expressed with these states the 
projection operator onto lead a reads 

Pa= [ de\e,a){e,a\. (12) 
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Using det6 = exp(trln6) in Eq. the logarithm of the characteristic function is given 
by 

ln(7(t,{A})=trln[l + 6(t,{A})]. (13) 

The thermodynamic hmit of the trace operation is defined in Eq. ( iTTl) . Mathematical 
subtleties concerning the existence of the determinant in Eq. ([8]) in the thermodynamic 
limit were discussed recently [10]. In the long time limit the "dot part" of the trace 
(see Eq. (11 II) ) gives a finite contribution which will be neglected in the following. 
Alternatively one can avoid the dot states alltogether by including them in (part of) the 
leads [8]. Using (e, a\nQ\e',a') = S{e — e')6aa'fa{^), where the /a(e) are the Fermi functions 
of the leads the matrix elements of the operator b{t, {A}) with the lead states are given 
by 

(ei,ai|6(t,{A})|e2,a2) = (14) 
Vrf(AJ / {e,,a,\e'^''j^e-'^''\e2,a2)dt'e-'^'^^fa,{e2) , 

where d{\a) = e*'^" — 1. Equations (fTSll and (fT4|) summarize the two tasks in order to 
obtain the characteristic function in the thermodynamic limit. In the first step the time 
evolution of the matrix element in Eq. f|T^ has to be determined. In the second step 
the trace in Eq. f|T3|) has to be performed. 

3. Long time limit 

As the current operators ja only involve operators localized near the dot region one can 
approximate the time evolution in the matrix elements in Eq. (fT^ in the long time 
limit by e~*^*'|e, a) e~*^*'|e, a+), where |e, a+) is the scattering state with outgoing 
boundary condition for the connected system [8] . Then the time integration in Eq. 
can be easily performed 

(ei,ai|6(t,{A})|e2,a2) ^ (15) 
2_]rf(Aa)(ei,ai + |ja|e2,a2+)e-^^"2/a2(e2) . 



Hei-e2j 

The fact that a simple analytical expression for the matrix elements has been derived 
allows systematic approximations for calculating the trace in Eq. (1131) . The first step is 
a formal power series expansion 

ln^7(t, {A}) = ^^tr 6"(t, {A}) . (16) 

n=l 

In the evaluation of the trace of 6" the time dependent prefactor of the sum in Eq. f|T5l) 
plays a central role. The product of the factors (e^'^J-'* — l)/ej,i with ej^i = ej — ei is used 
to obtain a product of n — 1 "energy conserving" delta functions 

gi.i,2t _ I gie2,3t _ 1 gi^n.it _ 1 ^ sin (ei,2t/2) sin (£2,3^/2) sin (e„,it/2) 

«ei,2 ^62,3 ' ' ' i^n,! 61,2/2 €2,3/2 ' ' ' en,l/2 

^ (27r)"-45(ei - e2)5(e2 - £3) . . . (5(e„_i - e„) . (17) 
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Therefore only one energy integration remains and the trace in the full one particle 
Hilbert space can be converted to a trace in the M- dimensional space of the lead indices 
a denoted by tT(^M) 

tr {A}) ^^Jde tr(M)c"(e, {A}) , (18) 

where the M x M matrix c(e, {A}) has the matrix elements 

Caia2(e,{A}) = 2nJ2d{Xa){e,ai + \ja\e,a2+)e-'^''^ faM- (19) 

a 

In contrast to the current matrix elements off-diagonal in energy the diagonal elements 
in Eq. (12 ip can be simply expressed in terms of the scattering matrix Sa^^aj 

(e). As 

shown in appendix A 

2n{e, ai + |j«|e, 02+) = sl^^{e)s^^^ (e) - 6a,aSa2a (20) 

holds[Ilj. If we define the M x M matrices e({A}) and /(e) as eaa'ii^}) = e*^"5aa' 
and faa'i^) = fa{^)Saa' the matrix c takes the form (suppressing the energy and A 
dependencies) 

c= {s^ese^ -l)f = cf . (21) 

Now the relation lndet(l + c) = trln(l + c) can be used backwards. With Eqs. (fT6ll 
and (JTSll this yields in leading time order 



lnc/(t,{A}) = — / de lndet(l + c) . (22) 

271 J 

Sub leading corrections increase only logarithmically with timepTl [HI [12]. The leading 
order term was correctly given by Levitov and LesovikfS] without presenting a derivation. 

4. Evaluation of the determinant 

In order to obtain explicit results for the leading order in time result for In g{t, {X}) the 
determinant in the integrand of Eq. f[22|) has to be calculated 

D{e, {A}) = det(l + cf) = det(l - / + s^ese^f) . (23) 

In both representations one has to calculate a determinant of a matrix which ist the 
sum of a diagonal matrix and an arbitrary matrix. The second decomposition in Eq. 
(!23l) is usually taken as the starting point [31 [13]. If one is interested in the probability 
distribution of the transferred charge in a single lead or the joint probability distribution 
of only a few (e.g. two) of the Qj it turns out to be preferable to use the first 
decomposition in Eq. (!23l) . Then we can use 

det(l + c) = 1 + tr c + ^ det(2)c(2) + J] det^^'^c^^) 

{i<j} {i<j<k} 

+ + det c , (24) 
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where e.g. det^^-'c*^^^ denotes a 3 x 3 sub determinant of c with the indices given by the 
summation variables. Because of c = c/ the Fermi functions can be factored out and it 
is sufficient to consider the subdeterminants of c. The matrix elements of c are given by 

W = - ^>ba' = Yl -ibdba'Sba' ■ (25) 

b b 

We begin with the discussion of the transferred charge in a single lead. 
4.1. Generalized Levitov-Lesovik formula 

We choose a = 1 as the channel index of the observed charge transfer, i.e. Ai is different 
from zero, but all Aq with a > 1 are put to zero. For a' > 1 this implies 

b 

while for a' = 1 one obtains 

Cai = ^ib^lsbi = {Sal - siiSu)dl . (27) 

b>l 

Apart from the additional term in cn the columns of the matrix c are all proportional 
to s^i- Therefore all subdeterminants of order larger than two on the rhs of Eq. ( !24l) 
vanish. Using 1 — = J2a=ii l-^iaP trace term reads 

tTc = Y\sia\\dlf\ + difa) . (28) 

The only non-vanishing 2x2 matrices are from the pairs {i,j) = (1, J > 1). In the first 
column the part proportional to s^j^ does not contribute and with didl = —{di + dl) one 
obtains 

det(2)c(i-) = -(di + dt)|si„|Vi/a- (29) 
Therefore the determinant appearing in Eq. (l22|) is given by 

det (1 + c) = 1 + 5^ \Sia\' [dlfafl + dlfja] = I + L, , (30) 

where fa = ^~fa- This is the generalized Levitov-Lesovik formula for the leading time 
contribution to In (^(t, Ai, 0, 0). The derivation presented here simplifies an earlier one 

4-2. Joint probability distribution for two leads 

In this subsection we present the general expression for arbitrary values of M for the 
characteristic function necessary to calculate the joint probability distribution for two 
observed charge transfers. We evaluate det(l + c) for arbitrary values of Ai and A2 and 
put all Aj for z > 2 to zero. For j > 2 this implies 

Cij = Y 4ii'^'^' - = slidiSij + sl2d2S2j ■ (31) 
I 
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The third and all higher columns of the matrix c are proportional to the two column 
vectors and s|2- Only the first two columns have an additional contribution. Using 
d2i — dl = d2 + d2d\ and di2 — d2 = di + did2 they are given by 

Cji = d*i5ii — s\id*iSii + s\2{d2 + d2d\)s2i , 

Ci2 = dl5i2 + s\^{di + didl)si2 - ^20^2^22 • (32) 

The structure of the matrix c implies that only determinants up to order four on the 
rhs of Eq. (12^ can be different from zero. 

Of the forth order determinants only those of the c^^^*-^-' with j > i > 2 are 
nonvanishing. As new objects 2x2 determinants formed by scattering matrix elements 
appear 

Sij^ = SiiS2j — SijS2i ■ (33) 

For j > i > 2 they contain the interference effect in the scattering of two fermions from 
the leads i and j to the leads 1 and 2. Using again d^d* = —{di + d*) one obtains 

^^^(4)^(12.,) ^ ^ ^ ^.^^gm\2fj^fj^ (34) 

The evaluation of all the determinants of order three and smaller is straightforward but 
tedious. It is therefore presented in appendix B. With the abbreviation 

B,^l{\S[f\'-\S^'\'-\s2,\' + M') (35) 
the general result for D{e, Ai, A2, 0, 0) = det(l + c) reads using the definition in Eq. 

m 

D = l + L, + L2 + [d,d;hf2i\s,2\' + Yl ^^•^^■) + (1 2)] 

jyi,2 

+ ^ E l^if + d*M){d2f2fj + d;f2fj) . (36) 

ij>2 

In the double sum the restriction to i ^ j is included by the fact that the S^^ vanish. 

The generalization to joint probability distributions of more than two charges is 
obvious but the expressions become rather lengthy. Therefore this will not be discussed 
further here. Instead we next present the result for all Aj 7^ for a special form of the 
scattering matrix. 

4-3. Separable scattering matrix 

The special case of at dot consisting of a single level, called "simple star" is described 
by a separable scattering matrix of the form (see Appendix C) 

Sij = 6ij + aiOjU . (37) 

The unitarity of the scattering matrix implies 

u + u* + \u\'^Y\ai\^ = . (38) 
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In the complex w-plane this is the equation of a circle with radius around (— rQ,,0), 
where l/r^ = Yli |ajp. For i ^ j this implies the inequality 



4| 




2 




to 






2' 


^1 



^^^■1 < 12 2 (39) 



used later. 

For the calculation of D(e, {A}) it is useful to write c as 

c= {shs - e)e^f = {shs - e)f . (40) 

Using the unitarity relation Eq.( !38l) and the definition 

S^ = J2\aiWe'^'-e'^^) (41) 
I 

the matrix elements of b defined by aba* = s'^es — e, i.e. factoring out the diagonal 
matrices a and a*, are given by 

bij = {e'^^ -e'^^)u+\u\^Sj . (42) 

It is easy to see that all determinants of submatrices of dimension three and larger 
vanish. For the m x m submatrices i and j take values ii with ii < i2 < ... < im- The 
determinant is unchanged if one subtracts the first row from all other ones 

-&n. = (e*^'-e^^n)«. (43) 

As the second and all higher rows are proportional to each other the sub determinants 
of dimension m > 3 vanish. The 2x2 determinants are readily calculated. Using the 
unitarity relation Eq. (!38|) the result for D takes the simple form 

D{e, {A}) = 1 + E \c^^\>J\M%f^f) = 1 + E - IW, 

= 1 + ^ + ^ \s,,\^d4]hf, (44) 

i i^j 

For the special case of the separable scattering matrix the interference terms discussed 
in subsection 4.2 and all higher ones vanish. 

5. Shot noise 

In this section we mainly elucidate our result Eq. (l36l) which determines the joint 
probability distribution for two leads and present separately results for M = 3 and 
M > 4 in the zero temperature limit. As we present explicit results we have to specify 
the Hamiltonian introduced in Eq. ([2]). In the one-particle Hilbert space the leads are 
described as nearest neighbor hopping chains 

hoa = -^i\m,a){m + l,a\ + }l.c.) . (45) 
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Two types of dots are considered. For the "simple star" the dot hamiltonian and the 
couphng term are given by 

M 

/i^ = K)|0)(0| , t; = -5^(rJl,a)(0| + H.c.) . (46) 

a=l 

For the dot consisting of a ring of Mdot sites pierced by a magnetic flux we have 

/.dot ^ J2 V,\z){z\ + J2{Ti^i+l\^){^ + 1| + H.c.) , (47) 

1=1 1=1 

where Mdot + 1 corresponds to 1. The coupling to the leads is assumed to be 

M 

v = -J2i^a\l,a){t,\+ll.c.) . (48) 

a=l 

5.1. Y-junctions 

For M = 3 the double sum in Eq. fl36p is missing and defined in Eq. 0351) takes a 
simpler form as |5'23'*P = |s3ip and |5'j3'*p = |s32p holds [3]. Using the unitarity of the 
scattering matrix Eq. (1361) simplifies to 

D = 1 + Li + L2 + [didlhh{\si2?h + \s2i?h) + (1 ^ 2)] . (49) 

For the three leg simple star |si2p = |s2ip holds and Eq. fHOj) reduces to the M = 3 
version of Eq. (jHl) with A3 = 0. At zero temperature Eq. ( H9l) reduces to D = I+L1+L2 
when lead 1 and 2 have equal chemical potentials. The correlation between the charge 
transfers to leads 1 and 2 is still present as hig is determined by \n.D. 

The joint probability distribution shown in Fig. 1 is for the case /^s = + A/x. 
If the energy dependence of the scattering matrix elements can be neglected in the 
energy interval [/i2,/^3] the energy integration in Eq. ( l22l) can be carried out and the 
characteristic function g is given by 

g(t,X^,\2,Q)=[l + \s^,\\e'^^-l) + \s2,\\e'^--l)Y' , (50) 

where Nt = tA/i/(27r). For integer values of Nt the A-integrations in Eq. can easily 
be performed and the Pnmit) in 

w{t,Qi,Q2) = ^Pnmit)S{Qi - n)5{Q2 - m) (51) 

n,m 

are given by 

Pn^{t)=( ^ V'' + "'V«r2™(l-Ti-r2)^-("+-), (52) 

\n + mJ \ m J 

where the Tj = |sj3p are the transmission probabilities from lead 3 to leads 1 and 2. 

For the "simple star" with Vq = 0, ri = r3 = 1 the transmission probabilities in 
the middle of the band are given by (see Appendix C) 

4 4|r2p 
^^=(2TN¥' {2+\r2\-y ^^^^ 
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For T2 = the connected system corresponds to an infinite ideal chain with perfect 
transmission Ti = 1. If the couphng to lead 2 is switched on Ti decreases and equals 
T2 = 4/9 for equal coupling T2 = I. 



t A n / (2 It) = 20 

- T3-1,X2-1/2 



O 10 



• • • • 

• • • • • 

• • • • • • 
• • • • • 



• • I 



10 
Q, 



Figure 1. Joint probability distribution for a three leg simple star for Nt = 20 
and three different coupling strengths T2 to lead 2 using the leading order time 
approximation for lug. The area of the dots is proportional to the p„.,„ in Eq. (j52p . 



Fig. 1 shows the probabilities Pnm{t) for Nt = 20 and three values of T2. For T2 = 
the open circle shows the "vanishing shot noise" which is an artefact of the leading 
order in t approximation for \ng. The neglected logarithmic corrections convert the 
single delta function to a Gaussian with a width proportional to (InA^^)^/^ for A^^ ^ 1 
fn] [HI [12] • For a weak coupling to lead 2 {t2 = 0.5, dark dots) the anticorrelation 
between the transported charges to leads 1 and 2 is clearly visible. It was discussed 
previously on the level of moments [Tlj. The light dots correspond to the case of equal 
couplings to the star. 

5.2. Systems with M > 4 leads 

For Y-junctions there are no interference effects in two particle scattering processes. 
The double sum term on the rhs of Eq. ([36]) only contributes for M > 4. To simplify 
the discussion we here discuss only the special case ni = ^2 = ^r, i.e. /i = /2 = fn at 
zero temperature which implies = Then Eq. ( 136|) simplifies to 

D = 1 + L, + L2 + ^Y1 \S'S\\did2fRhfj + dldlfnfjj) ■ (54) 

Already for M = 4 there are various possibilities for the chemical potentials of the leads 
3 and 4. The simplest one is to assume them to be equal. For arbitrary M > 4 and 
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f^i = fJ'L = I^R + ^IJ' for alH > 2 Eq. (IMI) further simplifies for the energy interval 
/iij < e < to 

D = l+Y,diY,\sij? + did2Y, (55) 

i=l,2 j>2 j>i>2 

The interference terms in the double sum also occur in the distribution of the total 
charge Qi + Q2 in the leads with chemical potential which can be obtained by 
putting Ai = A2 = A. As the prefactor d"^ = (e*'^ — 1)^ has no contribution linear in A the 
interference terms only enter the cumulants with i >2. When the energy dependence 
of the scattering matrix for /i/j < e < /xl is neglected the probability distribution for 
integer Nt = tA/i/(27r) ^ 1 is given by 

w'^'\t, Qi + Q2) = J2pi,'\t)6{Q, + Q2~n) (56) 

n>0 

with 

Pn^ = E L- /) / ) " ^ + 5)^-"+'5'(A - 2i?)«-2', (57) 

/ — Iq 

where Iq = max(0, n — Nt) and A and B are given by 

^ = ^:El^^.■^ B=j2\ss^r (58) 

1=1 j>2 j>i>2 
(S) 

In Fig. 2 we show results for the pn for a symmetric ring dot with = 0.5, Va = and 
Tj^j+i = e**/*^*^, where $ is the magnetic flux through the ring. We choose M = = 4 
with each ring site connected to one lead. The scattering matrix elements needed are 
presented in Appendix C . In addition to the probability distribution w^^^ of the sum 
Qi + Q2 (filled symbols) we show the probability distributions for the charge transfer 

to the single leads 1 and 2 (open symbols) which follows from Eq. fl30l) . The circles 
correspond to the case when the chemical potentials are close to the center of the band. 
Because |si2p = |si4p holds for e = for all values of $ the probability distributions for 
the charge transfer to leads 1 and 2 are identical. In the generic case they are different 
for $ 7^ 0(mod27r) as shown for e = — 1 (circles). For the parameters chosen it is clearly 
visible that the sum of the widths of these single charge probability distributions is 
smaller than the width of w^^\ This is another manifestation of the anticorrelation 
effect P mentioned in the discussion of Fig. 1. 

The joint probability distribution w(t, Qi, Q2) corresponding to the parameter val- 
ues used in Fig. 2 is shown in Fig. 3. In contrast to Fig. 1 where the asymmetry in 
the (dark dot) distribution resulted from asymmetric couplings, the asymmetry of the 
e = — 1 distribution in Fig. 3 is due to the magnetic flux. 

As a simple application of Eq. (jHj) for the separable scattering matrix we consider 
the distribution of the total charge Qr = Qi + ... + Qmr transferred to the first Mr 
channels which are all assumed to have the chemical potential /i/j. For the case where 
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^. tAn/(27i)=80, |Tii^;i=1 

it *=ll/10, E=0 ■ 




20 40 60 80 100 120 140 



n 

Figure 2. Probability distributions for a four site ring pierced by a magnetic flux 
$. The full symbols present the pn of the probability distribution of the sum 
Q — Qi + Q2- The open symbols show the probability distribution for the charge 
transfer to a single lead. The lines are a guide to the eye. Because of the magnetic 
flux the distributions for lead 1 and 2 (the two curves with open circles) differ except 
in the particle- hole symmetric case (open squares). 
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Figure 3. Joint probability distributions for a four site ring for the paramters used 
in Fig. 2 



the other = M — Mji channels have chemical potential /i^ = + AyU Eq. 0441) 
reduces to the M = 2 Levitov-Lesovik formula Eq. fl30l) with the effective transmission 
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probability 

Mr M 

^«ff = E E l^'^f<l- (59) 

i=l j=Mii+l 

The inequality follows from Eq. fl39|) or the explicit result for the simple star in appendix 
C. 

6. Summary 

A simple derivation of the leading time order result for the logarithm of the characteristic 
function which determines the full counting statistics for systems of noninteracting 
fermions was presented. The energy dependent determinant involving the scattering 
matrix of the M-lead system was simplified analytically for three cases. For the 
distribution of charge transfer to a single lead only the absolute values of the scattering 
matrix elements enter. For the joint probability distribution of the charges in two leads 
interference effects in the scattering of pairs of particles become important for M > 4. 
For a separable scattering matrix which describes a simple star like geometry the M xM 
determinant was evaluated for joint probability distributions for an arbitrary number 
of observed charges. All interference terms vanish and only the absolute values of the 
scattering matrix elements enter as in the generalized Levitov-Lesovik formula. 

Explicit results for the probability distributions were presented for the simple star 
and a ring pierced by a magnetic flux at zero temperature. Various manifestions of the 
anticorrelation effects in the charge transfer to the observed leads were shown. 

Finite temperature effects will be published in a forthcoming publication. 



Appendix A. Proof of Eq. (l20l) 

In this appendix we present a proof of the relation for current matrix elements with 
equal energies [H]. 

The scattering states |e, a+) introduced in section 3 obey the Lippmann-Schwinger 
equation [15] 

|e, a+) = |e, a) + go{e + iO)v\e, a+) , (A.l) 

where v ist the operator which describes the connection of the leads with the dot region 
and go{z) = {z — Iiq)'^ ist the unperturbed resolvent. 

As the projection operators Pa commute with ho the current matrix elements are 
given by 

+ ba|e,a2+) = ^(e,ai + \ [Pa,v]\e,a2+) . (A. 2) 

Using Eq. (fT2|) for Pa the first term of the commutator reads 

{e,ai + \Pav\e,a2+) = de {e,ai + \e ,a){€ ,a\v\e,a2+) ■ (A.3) 
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With help of the Lippmann-Schwinger equation the overlaps of the unperturbed states 
with the scattering states can be written as 

(e, ai + |e', a) = 5{e - e')Sa,a + ■ (A.4) 

This yields 

(e, ai + \Pav\e, 02+) = (e, a|f |e, a2+)(5aia (A. 5) 

,{e,ai + \v\e' , a) {e' , a\v\e, a2+) 
e-e' -iO ■ 

The matrix element of vPa can be calculated correspondingly. In the second term in 
Eq. (1A.5I) the —iO in the denominator is replaced by iO. This leads to a delta function 
for the commutator and the energy integration can be carried out. One obtains 

(e, ai + \ja\e, 03+) = i{e, ai + \v\e, a)5a^a - a\v\e, a2+)5a,a 

+ 2n{e, ai + |f |e, a) (e, a|f |e, 02+). (A. 6) 

With the definition of the scattering matrix Saa'(e) for potential scattering [15] 

w(e) = Saa' - 27rz(e, a\v\e, a'+) = 6aa' - 2mtaa'{e) (A. 7) 

the validity of Eq. (!20l) directly follows. 

The unitarity of the scattering matrix can be explicitely confirmed for the multi-lead 
system [8]. 




Appendix B. Derivation of Eq. ( I36p 

In this appendix we discuss the m x m subdeterminants of the matrix c defined in Eq. 
( I2TI) for Aj = for alH > 3. They are characterized by the index set in € [1,M] with 
ii < 12 < ... < im- The corresponding submatrices are denoted by c^*i'*2v,«m)_ ^j^g 
matrix elements of these submatrices can be directly read off the matrix elements of the 
full M X M matrix c 

Cij = Sii6ijdl + 6i262jdl + sl^aj + sl^Pj , (B.l) 

where the aj and jSj are given in Eqs. ( |3T1) and (|32l) . Like the full matrix c all columns 
of the submatrices are linear combinations of at most the four column vectors 6ii, 6i2, s]^ 
and sjg- Therefore all subdeterminants with m > 4 vanish. In the following we discuss 
the cases m = 4,3,2 separately, 
m = 4: 

Only for the index set {l,2,i,j) with 2 < z < j all four different column vectors 
appear. Therefore all other subdeterminants vanish. For the determinant of c^^'^'*'-'-' the 
following formula is used. Let a = and b = (■5J25 ■5j2)"^ two column vectors. 

Then the 2x2 determinant of linear combinations of the two vectors is given by 

det(aaj + bPi, aaj + bjSj) = {ai(3j — aj/5j)det(a, b) 

= {ai(3j - aj(3i){aibj - a^bi) . (B.2) 
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The use of this relation and Eqs. (PTi) and yield Eq. 
m = 3: 

All subdeterminants of d'^^ with 2 < i < j < k vanish. For i = 1 or i = 2 and 
2 < j < k the calculation proceeds using Eq. (]B.2|) 



detc^'^'^) = did2d;\sjiSk2 - SkiSj2\^ = -di{d2 + d*^)\sf^\'' ■ (B.3) 

In c'-^^*^ with i > 2 the first two columns Z = 1, 2 are of the type aia + jSih + (i;*e/ were 
is the /'th unit column vector while in the third column a^a + (3ib this additional 
contribution is missing. Therefore the determinants det(e/,a,6) and det(ei, 62, a(&)) 
have to be evaluated. This leads to 

detc(^2*) = (rfi + dl)dl{\S^^^\'' - + (1^2). (B.4) 

For m = 2 arbitrary combinations oi i < j in detc*^*-'^ contribute. The evaluation is 
straightforward using Eq. (lB.2p . 

In order to obtain the form presented in Eq. (!36l) the following "sum rules" for the 



were used which follow from the unitarity of the scattering matrix 



t>2 



j>i>2 

Appendix C. Results for the scattering matrices 

In this appendix the scattering matrices for the two dot models presented in section V 
are calculated. The t-matrix defined in Eq. (1A.7P can be expressed via the resolvent 
operator g{z) = {z — h)^^ as [I5] 



4a'(e) = (e, a|f |e, a') + (e, a\vg{t + iQ)v\t, a') . (C.l) 

For both dot models discussed in section 5 the first term on the rhs vanishes. With 
|(e,a|l,a)|2 = ^1 - (e/2)V7r this yields 

= 6aa' -2lTaTa'^Jl-{e/2Y{la\g{e + lQ)\la>). (C.2) 

For the simple dot all \ia) are given by |0) and the scattering matrix is of the separable 
form discussed in section 4.3. 

Using the projection onto the states on the ring the full resolvent matrix elements 
in Eq. flC.2p can be written as [S] 

{la\g{z)\la>) = {la\[z-K°' -l{z)]-'\la>) (C.3) 

with 

M 

l{^)=9l{z)Y,\ra?W){ia\ , (C.4) 

a=l 
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where /iq"* is the one particle Hamiltonian of the ring and g^{z) is the diagonal element 
of the resolvent of the semi-infinite chain at the boundary. For z = e ± zO it is given by 



g',{e±^0) = {eT^VA^)/2 . (C.5) 
The matrix inversion in Eq. (]C.3I) is trivial for the simple star because M^ot = 1 

For a symmetric ring with M = M^ot, Ta = t, Va = Vq, and Tj^j+i = TdotC**''^ the matrix 
inversion is straightforward using plane wave states on the dot [16] 



lrM^)\ln) = MY.^_y^_ |^|2^0(^) + 2rdot COS {k, + $/M) 



i=i 

where kj = 27ij/M. 
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